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Abstract: We propose a procedure to clarify part of the physical sector in the five di- 
mensional bubble geometries based on ideas similar to the split attractor flow conjecture 
OhI proposed by Denef. This procedure involves building some simple tree-like graphs that we 

r"! \ call skeletons without referring to the moduli space. The skeleton (tree) exists if and only 

if it passes the existence conditions which are purely based on some local CTC's (closed 
CN ' timelike curves) checking. Then, we propose the conjecture similar to Denef 's version which 

^. • states that every existing skeleton (tree) should correspond to some solution in which the 

global absence of CTC's is ensured. Furthermore, we propose two pictures to identify this 
correspondence explicitly and use some numerical examples to show how this procedure 
works. We also analyze the physical sector of the simplest bubbled supertube and see how 



^^ ' the existence conditions constrain the charge parameter space. 
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1 Introduction 

The black hole microstates counting [1] is considered to be one of the most important 
achievements of string theory. However, the counting is done in the limit of weak string 
coupling and it is unclear what these microstates become in a strong coupling regime when 
the gravitational back-reaction can no longer be ignored. On the other hand, the very 
existence of the black holes brings another deep puzzle: the information paradox. To 
address this issue, Samir Mathur and his collaborators propose an interesting idea, the so- 
called fuzzball proposal (see [2, 3] for a review). Roughly speaking, the proposal states that 
a classical black hole should be thought of as an effective description of the ensemble of many 
horizonless fuzzball states. In the stronger version of the proposal, the majority of fuzzball 
states can be realized as the solutions of supergravity. These solutions are completely 
smooth and horizonless and naturally can be considered as microstate geometries of a 
black hole. 

This proposal has achieved some great successes in the two-charge BPS solutions (for 
a review see [2]). However, to verify this proposal on the black holes with macroscopic 
horizons, one needs to extend the analysis to the three-charge solutions. There have been 
some developments in five dimensional BPS supergravity solutions [4-11] in this direction 
and a review is in [12]. It is clear that there are a large number of bubbled geometries which 
are smooth and horizonless. They are the natural candidates of microstate geometries of a 
given supersymmetric black hole or black ring. However, due to the difficulties of checking 
the conditions for the absence of closed time-like curves (CTC's), the physical sector of the 
bubbled geometries remains unclear. Without knowledge of this sector, it seems unlikely 
that a systematic study could be carried out or even to "count" the number of these 
microstate geometries. 

On the other hand, a similar problem is encountered in four dimensional multi-center 
BPS supergravity solutions but with the different physical origin. In [13, 14], an interesting 
criterion of the well-behaved solutions was proposed. This so-called split attractor flow con- 
jecture states that a multi-center solution exists if and only if there exist an corresponding 
attractor flow tree in the Calabi-Yau moduli space. There are also some studies of 4d-5d 
connection [6, 15-19] which suggest that the problem in four dimensional supergravity may 
relate to the CTC's checking in five dimensions. In [14, 20], it has been further suggested 
that this conjecture can be used in some five dimensional solutions as well by uplifting 
the four dimensional solutions. However, this uplifting procedure usually involves some 
rescaling of harmonic functions or coordinates and it is unclear whether there is always a 
unique way to connect four dimensional solutions to five dimensional solutions. Addition- 
ally, from the point of view of five dimensional theory, it is unclear why the solutions in five 
dimensional bubbled geometries should be related to the split attractor flow trees in four 
dimensional rather than five dimensional moduli space. And yet, the central charges and 
moduli space in five dimensional theory are all real and it seems unlikely to have the similar 
picture like the wall of marginal stability in the complex moduli space in four dimensional 
theory. 

To avoid the above subtleties and make the connection between the criterion in the 



conjecture and the solutions in bubble geometries more transparent, we take a more di- 
rect approach and propose two alternative forms of the conjecture based purely on CTC's 
conditions. These two forms are referred to as the "static" and the "dynamic" version 
respectively in this work. The "static" version is based on locating a "consistent" path in 
spacetime for a given solution while the "dynamic" version is related to some adiabatic 
property of the solutions of the bubble equations. Both of them are developed by consid- 
ering the "response" (i.e, decompose to two clusters) of the solutions under the tuning of 
the asymptotic constants similar to the original idea of the split attractor flow conjecture. 
However, to have a criterion that can be practically used to search CTC's free solutions, 
a systematic "way" of tuning the asymptotic constants must be specified . For that, we 
at first propose a systematic method to build the tree-like graphics. These graphics are 
referred to as the "skeletons" and are constructed by joining several single-center flows. 
Then, a criterion is proposed and it is amount to checking several CTC's conditions on 
these "skeletons". For those "skeletons" that pass the criterion, they will be one-to-one 
corresponding to some well-behaved components of the solution space of the bubble equa- 
tions. In fact, this criterion is very similar to Denef 's original split attractor flow conjecture 
even though it does not have the same physical picture related to the stability as we have 
explained. We will further explain the relation between the criterion, the two versions of 
the conjecture and Denef 's original conjecture in section 3.3. Although all of the different 
versions play some roles, it is the criterion based on the "skeletons" that will be the main 
tool in this paper. 

Unfortunately, the above picture based on the criterion does not cover all of the well- 
behaved components. Specifically, there are potentially large number of the components 
which simply resist breaking into two clusters even when the asymptotic constants hit the 
split point. These kinds of the solutions are referred to as the scaling solutions and for 
them, the picture will fail to distinguish different solutions or give an unambiguous answer 
as to whether the solutions are well-behaved or not. Nevertheless, through this split flow 
picture, one can explore the physical sector of the bubbled geometries at least partially or 
test the conjecture by checking the corresponding solutions explicitly. In this paper, we use 
the criterion to clarify the physical sector of the simplest bubbled supertube and also test 
the conjecture by examining several more complicated numerical examples. Additionally, 
we discover that in some systems there exists a bifurcation point between a scaling and 
a split flow component. In other words, the component of the solution space undergoes 
topological change and the split flow component emerges from the scaling component when 
the system is close to the split point. For more detail on this, see section 5.4. 

From a different perspective, this picture can be treated as some kind of coarse graining 
(or renormalization) scheme for solving bubble equations in which one reduces a large 
system of bubble equations to smaller and smaller subsets which become easier to solve. 
In other words, this basically collects some subset of the Gibbons-Hawing centers into 
a cluster and examines the CTC's conditions by treating this cluster as a single-center. 
Therefore, this might be precisely the scheme or positivity condition which are suspected 
to exist in [12]. However, as we have mentioned, the scaling solutions are like some kind 
of "entanglement" states formed by several centers and cannot be captured by this coarse- 



graining scheme. They are not understood well and may yet play very important roles 
in microstate geometries. Nevertheless, improving our understanding of the split flow 
solutions may provide the first step to resolve the mystery of the scaling solutions. 

The outline of this paper is the following. In section 2, we review the basic setup 
for bubbled geometries without fixing the asymptotic structures of the metrics. Then, we 
derive the bubble equations in general backgrounds. In section 3, we briefly review the 
wall-crossing and the split attractor flow conjecture. Then, we propose several alternative 
forms of the conjecture which are suitable to be used in bubbled geometries and explain 
the relation between them. In section 4, we explain the construction of the tree-like graphs 
( "skeleton" ) and their connection with the components of the solution space of the bubble 
equations. In section 5, we check several numerical examples. In section 6, we use the exis- 
tence conditions of the 3-center trees to identify the physical sector of bubbled supertubes 
and examining other exotic 3-center solutions. Finally, we present our conclusion and some 
discussion about potential future work in section 7. 

2 Bubbled geometries with unfixed asymptotic backgrounds 

2.1 Review of bubbled geometries 

It was shown in [4] that the general three-charge five dimensional BPS supergravity solu- 
tions can be obtained from the compactiflcation of the eleven dimensional M-theory. For 
T^ compactiflcation, the metric has the following form: 

ds?i = - (j^^) (dt + kf + {ZiZ2Z^f'^ h,nndx^dx^ + (^] {dxl + dx^) 

^) (dx'r + dxj) + (^) (dxl + dxU (2.1) 

A = A^ A dx5 A dxQ + A"^ A dxj A dxs + A^ A dxg A dxio , (2.2) 

where A and k are one-forms in the flve-dimensional spacetime. The base space metric 
hmn is a four-dimensional hyper-Kahler metric. The warp factors, Zj, can be obtained 
from solving three linear BPS equations sequently: 

v'^Zi = ^CijK M (e-^ A e^) 

dk + Mdk = ZjQ^ , (2.3) 

where the self-dual "dipole held strengths" @ are deflned as: 

e^ ^dA^ + d (^) , (2.4) 

where *4 is the Hodge dual taken with respect to the base space metric, hmn- 

To construct an explicit solution, one needs to specify a hyper-Kahler metric for the 
base. In [5, 6], the Gibbons-Hawking space was considered to be the base and it was shown 



that many completely smooth five dimensional solutions could be explicitly constructed. 
These are the so-called bubbled geometries. The Gibbons-Hawking space is a subclass of 
hyper-Kahler space for which the metric is explicitly known: 

h^ndx'^dx'^ = V{dx\ + dxl + dxl) + ^{d^) + A- dxf , (2.5) 

where V is an arbitrary harmonic function in M'^ and A is^: 

Vx A = VV . (2.6) 

The solution of BPS equations (2.3) in this base space is specified by eight arbitrary 
harmonic functions {V,K^,K'^,K'^,Li,L2,L3,M). The warp factors Zj are: 

1 K-^K^ 

Zi = -CijK^^^ + Li, (2.7) 

and the one-from, k, is: 

k = ^iidil^ + A) + oj, (2.8) 

where n is: 

and LO satisfies 

Vxo; = VVM - MVV + -{K^VLi - LiVK^) . (2.10) 

For T^ compactification, Cjjk = \^ijk\- To have smooth geometries, we need to further 
overlap the sources of these harmonic functions. In other words, for a system with N- 
centers, we have the following harmonic functions: 

N N U 

Li = {lo)i +Et^' M = ^0 + Y^, (2.11) 

where rj is the distance to the center 'j'. The regularity of Zj and fi at rj = implies {lj)i 
and rrij should be: 

1 7 7 -L 7 7 7 

\h)i = - -^CiJK , mj = —CijK 2 • (2.12) 

2 Qj 12 g| 

The constants of the harmonic functions, {q^jk^, (Zo)/,?7io)) will determine the asymptotic 
structure of the metric. We will leave them unfixed at this moment. 



^One should not confuse A with the three gauge fields, A^ , which are coming from the reduction of the 
three-form fields, A 



2.2 Bubble equations in general asymptotic background 

With the asymptotic metric unfixed, the parameters of the bubbled geometries with N cen- 
ters will include the eight asymptotic constants, the charge parameters (i.e., {qj ,k^ ,k'j , /c?)) 
and the positions of all centers. However, not every solution with arbitrary parameters will 
be physical. In general, these solutions will contain CTC's (closed-time-like curves) and 
require one to constrain these parameters carefully to avoid their appearance. Specifically, 
it was shown that the necessary conditions to be free of CTC's are: 

Q > 0, VZi > 0, 7=1,2,3, (2.13) 

where, 

Q = Z1Z2Z3V - fi^V^. (2.14) 

These conditions are the main sources of difficulty to understanding the physical sector of 
the bubbled geometries because one needs to check these conditions globally in order to 
make sure the solutions are well-behaved. However, there are some equations that can be 
obtained by verifying the first condition in (2.13) at the positions of the centers. These 
so-called bubble equations are relatively easier to check and are proven to be quite useful 
even though they alone cannot guarantee the well-behavedness of the solutions. 

To derive the bubble equation, recall that Zi are regular while V diverges at the centers 
and therefore it is clear that fi must vanish at the position of every center in order for the 
first condition in (2.13) to be satisfied. Using the definition of fx in (2.9) and harmonic 
functions in (2.11), the constant part of /i at each center is: 

jLl J,2 E3 _L h.2 7,3 El I l3 ul h,2 l1 1,2 r,3 

^{r. = = -^-^ ^-^ ^-^ 2 ^ i ^ q + m 

qj qj 

- ^ E (^ Hih)l - H(l)l - i^lk') , (2.15) 

zqj ^^^ \qj / 

where {q,k^ ,(l)j,fh) are constants including the contributions from the asymptotic con- 
stants and other centers. By using the regularity constraints in (2.12), after some algebra, 
one can show that the conditions, fi{rj = 0) = 0, can be concisely written as: 

-E^-^T^ = (r,'/i), i = l,2,...iV, (2.16) 

where N is the number of centers , rij is the distance between center 'i' and 'j' and Tj is a 
charge vector^: 

r, = [q,, k',^CjjK^, ^CjjK^^) . (2.17) 

Similar to charge vectors, h is an eight-component vector formed by the asymptotic con- 
stants: 

h = (go, A:^, -(/o)/, 2mo), (2.18) 



In bubble geometries, the charge parameters of the center can be specified by its magnetic part alone, 



for example, r'*" = (qj, k' 



and generally we will refer to it by the asymptotic vector. Finally, the symplectic intersec- 
tion product is defined as: 

{A\ A', Au A^),[B\ B\ Bi, Bo)) := - A^ + A' Bj - AiB' + A^B'' . (2.19) 



Notice that for the bubbled geometries, this symplectic product of the two charges is 
actually equivalent to the flux between them as defined in [5]: 

(r,,r,) = ^,<?,n«n;fng), (2.20) 

with n}- = ( ^ i . Additionally, for convenience, we define the vector, 7^, which is 

composed of the eight harmonic functions as the following: 

n = (v, K^ -Li, 2m) . (2.21) 



To fix the geometry to be asymptotically fiat, h must be chosen appropriately. If we 
normalize the constants of Zi to be one, h should be set to: 

3 N 

ho = (0, 0, 0, 0, -1, -1, -1,-Y.Y. ^i) • (2.22) 

1=1 j 

By fixing h at /iq and using (2.20), the bubble equations in (2.16) reduce to the familiar 
form as defined in [5]. However, one notices that for general h, the summation of all A^ 
equations will not be satisfied because the left-hand side of (2.16) sums to zero due to 
the antisymmetric property of the symplectic product (2.19) while the right-hand side is in 
general not zero. Therefore the self-consistency requires h to satisfy the following condition: 

N 

rt = ^r,,/i) = 0. (2.23) 

j 

This consistent condition of h and the general bubble equations (2.16) will play an impor- 
tant role in section 3. 

2.3 Large gauge transformation 

It was noticed in [17] that there is a symmetry of the solutions associated with the shifting 
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c^c-^K^). (2.24) 

This transformation does not change the metric, CTC's conditions or bubble equations but 
shifts the gauge fields by some constant forms: 

A^ ^ A^ + c^ {di/j + A) . (2.25) 



Since it changes the boundary value of the gauge field, the symmetry associated with this 
transformation can be regarded as some kind of large gauge symmetry. 

For a multi-center solution, this transformation is equivalent to shift every charge 
parameter and asymptotic constant as the following: 



kj -> kj + c^ Qj , 



qo,kl + c^ qo , -{lo)i + CijK c-^ k^ + - Cjjk c^ c^ 



2 



qo 



2 mo - c^ (/o)/ + ^ CijK {qo c^ c"^ c^ + 3 c^ c^ k^) ) . (2.26) 

Because all centers in bubble geometries have nonzero Gibbons-Hawking (GH) charges, 
gj, every ki will be shifted by some constant under this gauge transformation. Therefore, 
for every center, there exists a particular gauge choice such that the center contains only 
the GH charge (i.e., T = (g, 0,0, 0,0, 0,0, 0)). This implies that ki may not be considered 
as physical quantities since they are not gauge invariant. However, some combinations'^ 
or relations of ki for several centers may become gauge invariant. Therefore, to define 
the physical sector, one needs to have some constraints on this combinations or relations 
between centers rather than the exact values of ki . This is probably one of the reasons why 
it is difficult to define the physical sector in bubbled geometries. 

3 Split (attractor) flow conjecture 

We have explained the difficulty of checking CTC's condition globally in five dimensional 
bubbled geometries and it turns out that a similar difficulty is shared by four dimensional 
BPS multi-center solutions. In four dimensions, the metric is determined by the entropic 
function, S, and it can potentially become imaginary somewhere in the space. To ensure the 
well-behavedness of the solutions, one needs to check if S is real globally which is similarly 
difficult to the CTC's checking. To overcome this difficulty, an interesting conjecture was 
proposed in [13, 14]. The conjecture states that a multi-center solution exists if and only 
if a split attractor flow tree exists in the moduli space. A split attractor flow tree is a 
tree-like graph in the moduli space and it starts from the moduli at infinity and terminates 
at attractor points of each of the centers. The flow can only split when it crosses the wall 
of marginal stability which is some codimensional one surface in moduli space. On the 
wall, the phases of the central charges of two constituents are aligned and can be separated 
infinitely far away. 

In order to use this conjecture to explore the physical sector of the bubbled geometries 
or verify the conjecture, we need to identify the split attractor flow trees in five dimensional 
moduli space and find the explicit connection between the trees and the solutions. How- 
ever, the central charges and the moduli space in five dimensional theory are not complex 
and it seems unlikey to have a similar picture as with the wall of marginal stability in 
four dimensional moduli space. On the other hand, it was suggested in [14, 20] that one 
can study five dimensional solutions in a similar manner through the uplifting of 4d to 5d. 



Typical examples are the total asymptotic charges and angular momenta of the system. 



However, this uplifting procedure usually involves the rescaling of the coordinates or har- 
monic functions which seems to make the connection to five dimensional solutions become 
ambiguous. Also, the lacking of the understanding of how this conjecture works from the 
point view of the five dimensional theory is unsatisfactory. Here, we will take a more direct 
approach. We identify some tree-like structures from the five dimensional solutions without 
referring to the moduli space and propose the criterion of the existence of the trees and 
solutions solely based on CTC's conditions. We will still refer to these tree-like structures 
as the split flow trees due to their similarity with the trees in the four dimensional moduli 
space. We propose two methods^ to construct the pictures of the split flow trees. Then 
from these two pictures, we propose two forms of the split flow conjecture. Finally, we will 
compare these two versions of the conjecture with Denef 's original conjecture. 

3.1 "Static" split flow trees in space and the conjecture 

In this section, we present a method to "locate" the flows in spacetime once a solution 
is given. For bubbled geometries, the space where the centers are located is M^ inside 
the Gibbon-Hawking space. To see the image of the flow in this space, one just needs 
a simple observation that when changing the asymptotic moduli to any point along the 
flows, the system should become another valid solution with a different asymptotic metric. 
For a single-center solution with spherical symmetry, the flow is a single line in the moduli 
space. Any point on this line map to a sphere in M.^. Then, taking the asymptotic moduli 
to any point down the flow will result in another valid solution with a truncated flow. 
However, for a two-center solution, the situation is different because we don't have spherical 
symmetry but only an axial symmetry. So, what should be the trace the flow follows in M^? 
Specifically what is the orientation it should take when flowing? To answer this question, 
one notices that for an arbitrary path, the values of the harmonic functions at the points 
along the path will not satisfy the consistent condition (2.23). This means if we assign the 
asymptotic constants by the values of the harmonic functions at these points, there will be 
no solution for the bubble equations. This, of course, implies that the "image" of the flows 
in space cannot take any arbitrary path presuming that assigning the asymptotic moduli to 
any point down the flow in moduli space should result in another valid solution. Therefore, 
the correct path should be confined by the consistent condition. For a two-center solution, 
we have the following condition: 

(r,,?^(r)) = /rt,^ + ^ + /ioo) = {Ta,n){---) = 0, (3.1) 

where /iqo is the asymptotic vector and we have used the fact that (Ft, hao) = 0. One can 
see that if (Fq, F{,) ^ 0, this condition defines the equal-distance plane between two centers. 
Therefore, the main flow in the moduli space maps to a plane in M^ while a point on the 
flow maps to a circle on this plane due to the axial symmetry. Similarly, to find the flows 



''We use "static" and "dynamic" to distinguish two pictures. Sometimes, we ignore the quotation marks 
but it should be understood that they should not be interpreted literally. 
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Figure 1. This picture represents the two-center split flow in space and was obtained by cutting a 
slice along the axis of the two charges at a particular azimuthal angle. The solid blue and red lines 
represent the actual flows while the whole lines included the dashed ones are defined by (3.1-3.3). 



which terminate at Ta and r^, we need to consider the following conditions: 









+ (Fa, /loo) = , 

+ (Fb, /loo) = . 



(3.2) 
(3.3) 



One may notice how they are very similar to the two-center bubble equations. In fact, to 
solve these equations is equivalent to solving the bubble equation and it gives r^ = r^ = rab- 
Therefore, these conditions define two spheres centering at Ta and Ff, and with the radius 
equal to the separation of the two centers. To see this more clearly, we can look at a slice 
at a particular azimuthal angle and the flows will be look like Figure 1. The blue line 
represents the main flow while the red lines represent the split flows. If one parametrizes 
the blue line by the distance to one of the charges (i.e., r = Va = rb), then the domain of 
this line is r € [oo, -^J. The end of the line is at r = ^ and we will call it the "half-split 
point". The point where the blue and red lines cross (i.e., r = Vah) is called the split point. 

For the solution with just two centers, the tree-like pattern is clear. However, if we 
include more centers, we need to consider many possible lines with different combinations 
of charges. The resulting pattern may be quite confusing. The best way to extract the 
tree from the pattern is probably tracing the lines starting from the centers. When the 
two lines cross, remove the remaining lines and join them to the line corresponding to the 
summation of the two charges. Eventually, one trace it back to the line that goes to infinity 
and the result is the tree that corresponds to this solution. 
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Figure 2. These are several simple examples of the configurations which will violate the split flow 
conjecture if they exist. The yellow circle represents the region full of CTC's. 

With this picture of the spht flows in R^ , we can state the static version of the spht flow 
conjecture: for a solution in which a clear split-flow pattern can be identified by the method 
described above, the absence of the CTC's along the paths of the flows guarantees that the 
whole space is CTC's free. In other words, this version of split flow conjecture implies that 
the situations in Figure 2 cannot exist. Among them, the first one is particularly interesting 
because we know for some multi-center charge configuration, the attractor point of the total 
charge is not well-behaved. This usually means there will be CTC's near this attractor 
point and since the half-split point is closer to the attractor point than the split point, it 
seems possible that the half-split point dips into CTC's region while the split point remains 
free of CTC's. Therefore, testing CTC's near the half-split points for those solutions with 
clear split-flow patterns provides an interesting way to test the conjecture. 

This "static" version and the corresponding conjecture provide the representation of 
the trees in spacetime and a natural resolution if one wants to find out the tree that 
a particular solution corresponds to. However, there exists several drawbacks for this 
approach. First, from several numerical examples, we found out that if the fluxes that 
are involved in the split points of the tree are very different from each other in scale or 
very close to zero, some lines may be broken even though the solution itself is indeed a 
split flow solution based on the "dynamic" analysis in the next subsection. This not only 
makes it difficult to extract the tree from the solution unambiguously but also makes it 
hard to distinguish this solution from the scaling solutions in such cases. Furthermore, in 
order to use the "static" version of the conjecture, one needs to construct solutions first. 
Therefore, it is clear that this version is not suitable to be used in any systematic search. 
Due to these drawbacks, we will mostly rely on the "dynamic" version combined with the 
"skeleton" construction which will be introduced in sections 3.2 and 4.2. However, due to 
its close connection with spacetime, the "static" version is useful for testing the conjecture 
directly. Therefore, we will mostly use it to verify the conjecture by looking for any counter 
example. Specifically, we will check the CTC's conditions at the half-split points when we 
verify the existence of three-center trees in section 6. 

3.2 "Dynamic" version of the split flow trees and the conjecture 

Corresponding to the wall-crossing in the moduli space, what one will see in the four di- 
mensional supergravity picture is that when the asymptotic moduli (vector) approaches the 
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wall of marginal stability, the system will be divided into two parts which are separated 
further from each other under the constraint of the integr ability conditions. Since the bub- 
ble equations (2.16) in bubbled geometries have exactly the same form as the integrability 
conditions in four dimensional supergravity, it seems quite natural to follow this thinking 
and consider how a solution deforms under the constraint of the bubble equations when 
we change the asymptotic vector, h. 

To be more specific, consider a certain number of centers with the total charge equal to 
Ft. Then we start with some asymptotic vector, h^o, and begin to change it "continuously" 
by following the flow: 

hf = hoo + r]Tt, r/iO-^oo. (3.4) 

Then we "continuosly" trace how this solution changes under the constraint of the bubble 
equations. At some point, one may find the solution is divided into two parts which are 
separated further away when r] is approaching a particular value, 77. After t] passes this 
value, the solution no longer exists because it decays to two parts separated infinitely far 
away. However, one can trace the two parts separately starting from 77 and go on until 
eventually the whole solution reduces to several independent single centered solutions. If 
one find a solution can be completely deconstructed in this way, one can conclude it is 
a split flow solution and the order of the splitting determines the tree it corresponds to. 
However, in some cases, one may find some or all of the centers form some composite 
element which does not expand or separate to two parts no matter how one changes the 
asymptotic vector and on the contrary, it usually decreases in size when one follows the flow 
in (3.4). This kind of behavior is the signature of the scaling solutions and unfortunately, 
we cannot use the split flow techniques to search them or classify them. 

Notice that the solution is not uniquely determined by the bubble equations. What 
we trace here is just a point in one of the components of the solution space which will be 
explained more in section 4. However, this behavior of splitting should not be considered as 
the feature or property of a particular point in the solution space but the whole component. 
In other words, it does not matter which point in the component you are tracing^ as long 
as you do it continuously. The reason we should do it "continuously" is to prevent the 
solution jumping to the potentially many other components nearby. This indeed causes 
some difficulty in implementing this picture on the solutions practically. However, in 
general, it seems to provide a more accurate and unambiguous classification of the solutions 
compared to the "static" version. 

From this picture, we propose another form of the split flow conjecture as the follow- 
ing: The fact that if a solution contains CTC's or not cannot be changed by any continuous 
deformation under the constraint of the bubble equations. This includes those deformations 
induced by tuning the asymptotic vector continuously as long as one keeps it inside of the 
physical region. More precisely, what we mean is that the asymptotic vector is free to move 



^There are some subtleties if the scaling and split flow solutions coexist. In this situation, it is possible 
to have a bifurcation point between scaling and split flow components. If such a point exists, it will be 
crucial to decide which part of the component to trace. For more detail, see the numerical example in 
section 5.4 
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around or even split when it hits the spht point but it must be bounded inside of the re- 
gion where the corresponding self-consistent condition in (2.23) is satisfied and the metric 
at infinity determined from it must be CTC's free. While we have been considering this 
"dynamic" procedure involving tuning the asymptotic vector as some kind of operational 
definition of saying one solution belongs to a particular tree, it may have a physical inter- 
pretation as well. The nature guess is that it corresponds to changing the background of 
the whole system adiabatically. For example, the Taub-Nut interpolation in [17] that con- 
nects five dimensional black rings to four dimensional black holes can be regarded as some 
special fiow in this "dynamic" picture. Prom this point of view, the "dynamic" version of 
the conjecture seems quite natural. Because it is very hard to believe that one can render^ 
a well-defined solution unphysical by changing the background adiabatically presuming, of 
course, the background itself is physical. 

Notice that in the "dynamic" version of the conjecture, how one should tune the 
asymptotic vector was never specified. However, in order to use the conjecture practically, 
one needs to construct a systematic method to tune the asymptotic vector. In the beginning 
of this subsection, we have introduced a simple choice like the one in (3.4) which is a single 
centered fiow. In fact, one can construct a complete split flow tree solely from joining 
these single centered flows. The result is a tree-like graph referred to as the "skeleton" and 
its construction will be explained more in section 4.2. Moreover, we will also introduce 
a criterion for identifying the well-behaved trees based on checking CTC's conditions on 
these skeletons. Unlike the "static" version of the conjecture, we do not need to solve the 
bubble equations in order to construct a "skeleton" and therefore, this criterion is useful 
for a systematic search. 

3.3 The relation between different versions of the conjecture 

We have introduced several versions of the conjecture and here we clarify the relation 
between them and Denef's original split attractor flow conjecture. 

As we have mentioned in the previous subsection, the "dynamic" version never specify 
a way to tune the asymptotic vector. From the perspective of the "dynamic" version, to 
identify a split flow solution is equivalent to ask if there exist a way or path to tune the 
asymptotic vector such that the solution can be completely deconstructed. However, this 
path by no means needs to be unique. For example, in the "static" version, once the split- 
flow pattern is clearly identified and no CTC is discovered on the pattern, it immediately 
provide a valid choice of the path for the "dynamic" version since the harmonic functions 
are analytic and they satisfy the constraint in (2.23). Furthermore, the precise path of 
this pattern will be different depending on the internal moduli of the solution and all of 
them are valid choices of the path. Additionally, a "skeleton" that pass the criterion in 
the next section is also a valid choice. Therefore, it is clear that the "dynamic" version is 
stronger than any criterion based on a specific choice of path including the "static" version 
and the "skeleton". On the other hand, it is also clear that the "dynamic" version is more 



^Naturally, when we take the asymptotic vector to the spht point, the system will be divided into two 
parts. However, the resulting two independent systems are still both physically well-defined. 
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general because it is supposed to be applicable to any solution including a scaling solution. 
However, in scaling solutions, there may be multiple components including good or bad ones 
connected to the same scaling point at which all of the centers are overlapping with each 
other. Clearly, if the scaling point exists, this will violate the "dynamic" version because 
through this point one can connect a good component to a bad one. However, from the 
speculations in [9] and the detailed analysis in [21], it seems quite likely that the scaling 
point will be removed by some quantum effects. If this is true, the "dynamic" version 
can still be true even for scaling solutions. Even so, the "dynamic" version is not exactly 
a criterion but more like some kind of adiabatic property of the solutions. Therefore, It 
cannot help us sort out the scaling components because we cannot "disentangle" the scaling 
solutions by tuning the asymptotic vector. 

Even though there exist several differences between the different versions, it is still 
possible that they are all equivalent when all of them are applicable. One can image that 
the several choices of path including those in the "static" version and the "skeleton" form 
a representative "bundle" of paths. That is if any path in this bundle hit any pathology, 
there is no other choice of path can do better job. If this is true, all of the different versions 
will give same verdict to any solution when all of them are applicable. This is the point 
of view we take in this paper. On the other hand, it is also possible that the "dynamic" 
version is too strong and only the criterion based on the certain path choice is true. In that 
case, the difference between the "static" version and the criterion based on the "skeletons" 
will be an interesting open question. 

If we accept the above "bundle" picture, it is quite natural to select the "skeleton" 
as the representation of the whole "bundle" because it can be systematically constructed. 
Therefore, the criterion based on the "skeleton" can be used to search the healthy "bundles" 
and the corresponding trees. As we have mentioned, the "skeletons" are constructed by 
joining several single centered flows. If we change the setting to four dimensional theory 
and map the asymptotic vector to a point in moduli space, a "skeleton" will look like a 
split attractor flow tree. If we change the physical conditions accordingly, for example, 
instead of checking CTC's conditions, we make sure the central charge never hit zero in 
the middle of the flows, etc., the criterion in the next section will be exactly like Denef's 
original conjecture. This suggests that the criterion should be considered as an equivalent^ 
or a natural extension of Denef's original conjecture to five dimensional theory. 

4 Split flow trees <^=^ split flow components 

4.1 Solution space of bubble equations 

Since a solution must satisfy the bubble equations in order to be a candidate for well- 
behaved solutions, it is essential to understand the solutions of bubble equations. To 
be more precise, we should call them a solution space. The reason is that the bubble 
equations cannot fix the solutions completely because the number of equations is smaller 



^However, this will imply the causality (e.g. CTC's conditions) in five dimensions are connected to the 
stability issues or the wall-crossing phenomena in four dimensions. How that is so may be an intertesting 
question worth to explore. 
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than the degrees of freedom from the positions of the centers. Therefore, the solutions have 
some moduU or the "shding" degrees of freedom left to deform themselves. In [21], the 
space spanned by these moduli was referred to as the solution space^ and was studied in 
some simple cases. In the paper, the symplectic structure of the space was identified and 
the solution space was promoted to the phase space and quantized. By doing so, the wall- 
crossing formula [14] was recovered. This suggested that the wall-crossing or split attractor 
flows may have a very close relation with the solution space of bubble equations. However, 
it is still difficult to say something concrete about this space especially when the effective 
number of pairs with nontrivial fluxes between them is not so small. Nevertheless, through 
numerical studies, it is clear that in general the solution space is not even connected. 
This means that in most situations, the solution space has multiple components and each 
component is spanned by its own "sliding" degrees of freedom. Here, we roughly classify 
the components of solution space into three categories: 

• Split flow components: Any solution in these components can be deconstructed 
to several single-center solutions by tuning the asymptotic vector continuously along 
the split flow trees. There exists a one-to-one correspondence between these kinds 
of components and split flow trees. They are the main subjects of the split flow 
conjecture and their absence of CTC's globally is guaranteed by the conjecture. 

• Scaling components: the solution of this kind has some or all of its centers forming 
scaling solutions. Roughly speaking, these centers trap in some kind of "entangle- 
ment" state and cannot be torn apart by tuning the asymptotic vector. The whole 
solution may still have tree-like structure with some flows terminating at the scaling 
solutions. From the perspective of the split flows, these scaling solutions should be 
treated as single-center points or attractor points and should only exist if the cor- 
responding attractor points are well-behaved. Because there are potentially many 
scaling components corresponding to the same attractor point in the bubbled geome- 
tries, the one-to-one correspondence is no longer true in these cases. Furthermore, 
if a solution shows scaling behavior, the conjecture can no longer help us determine 
whether this solution is well-behaved or not. Depending on whether the solution is 
truly absent of CTC's globally or not, they will belong to the scaling components or 
the bad components. 

• Bad components: the solutions of this kind contain CTC's even though they satisfy 
the bubble equations and should be regard as unphysical. 

From the perspective of microstate geometries, the scaling components are quite interesting 
and probably very important. Unfortunately, even though there are some numerical exam- 
ples showing their existence, there is no systematic method to search or even distinguish 
the scaling from the bad components. On the other hand, the split flow solutions which 
belong to the split flow components are under better controls and we will show how to find 
the explicit connection between them and split flow trees. 



*In this paper, it is tlie solution space of the integrabihty conditions that was studied but the conditions 
are practically the same as the bubble equations in (2.16). 
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4.2 Building the "skeletons" of the split flow solutions 

In this section, we demonstrate the procedure to identify all split flow components given 
a set of centers with charges, Fj, and the asymptotic vector, /iqo- At first, we exhaust 
every possible split flow tree and build the "skeleton" for each of them. A "skeleton" 
is a tree-like graph that contains the complete information of all split points, how the 
harmonic functions change along the flows, and finally, the topological data of the tree. In 
other words, a skeleton can be seen as the blueprint of a solution. Then, we propose the 
existence conditions for the skeletons and use them as a criterion to pick out the physical 
trees. 

As an example, we will use a simple three-center system and will adopt the convention 
in [14] to use a nested list for identifying a tree. To start, consider a set of three centers, 
{Fa , Tb , Tc}, there are three possible trees: 

{r,,{r;,,r,}}, {rfe,{r,,rj}, {r,, {r„, r^}}. (4.i) 

From a given asymptotic vector, hoo, the structure of the tree and the charge of each center, 
we can construct the "skeleton" for each tree. For example, the skeleton of the first tree 
is shown in Figure 3. One can see that every edge in the skeleton looks like an attractor 
flow of a single-center solution and every split is treated like a two-center split. How the 
harmonic functions change along the flow and the positions of the split points are therefore 
determined accordingly. Every flow is parametrized by some positive variable, r/j, which 
is related with the inverse of the radius of the corresponding single-center solution. The 
parameter, r/j, starts from zero to infinity if the flow terminates at a center or stops at some 
particular positive value, r}j, if it hits a split point. The exact value of f/,, for a particular 
split is determined by the flux between two centers and the value of the asymptotic vector 
at the beginning of the flow. For example, fjt and r/^c are: 

m = -^ = _ _iIiLli^22l^ , (4.2) 

^^c - ^ - -^f^TTJ- (4-3) 

It is immediately clear that not every skeleton exists because rji must be positive while it 
is not necessarily true for arbitrary charges and h. Also, in order to use the skeleton as the 
path for tuning the asymptotic vector just like the "dynamic" procedure in section 3.2, we 
need to ensure there is no CTC along every edge in the skeleton. Therefore, we have the 
following existence conditions: 

• The positivity of the separations: rji > for every split point. 

• Absence of CTC's along the flows: Q (r/i) > , VZj {j]i) > . 
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Figure 3. This figure shows the "skeleton" of the tree, {Ta , {Ti, , Fc}}. The edges are like the singe- 
center flows parametrized by rj. How the harmonic functions change along the flows is completely 
determined once the asymptotic vector, hoc, is chosen. The two split points are present and labeled 



IS hsi and hg^ ■ 



The conditions above technically reduce the CTC's checking from three dimensional space 
to several one-dimensional lines. However, it is still difficult to check the second condition 
in practical and we adopt a weaker form of the existence conditions which reduce the 
checking to several points: 

• ?7i > for every split point^. 

• Absence of CTC's at every vertex and end point, namely: 

Q(r,)>0, VZi{Ti)>Q, Qihi)>0, VZi{hi)>0. 

Here we would like to explain more about Q{Ti) > and VZf(Ti) > 0. In bubbled 
geometries, for a true single-center flow that starts from some h and terminates at a GH 
center, the CTC's conditions are greatly simplified. For example, for a flow starting from 
hs and terminating at F'"^^ = (q, k ), the harmonic functions change according to the 
following: 

n = Trj + hs. (4.4) 



^One finds that this condition is related to the existence condition of the wall of marginal stability in 
[14]. The subtle difference arises when the wall of anti-marginal stability is involved. However, it is possible 
that the anti-marginal stability will be ruled out if one includes the second condition. 
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As r is a single GH center, we can always choose the gauge such that r^™' = (g', 0,0,0). 
Supposedly, /is becomes h'^ in this gauge: 

K = iqo,ki,-ilo)i,0), (4.5) 

where rriQ is required to be zero by the self consistent condition (i.e., (F^™', h'^) = 0). Then 
VZj and Q are greatly simplified: 

VZi{7]) = qilo)ji] + 0{1), 1=1,2,3, (4.6) 

Qiv) = q{lo)i{lo)2{lo)3V + 0{1). (4.7) 

Since all of these functions change linearly with rj, we only need to check the beginning 
and the end to ensure the whole flow is free of CTC's. The starting point is either at 
infinity or some split point and should already be verified in the other existence conditions. 
The conditions at the end can be imposed simply by requiring that the slopes in (4.6) are 
positive: 

(?(/o)i>0, q{lo)2>0, q{loh>0. (4.8) 

Although these conditions are expressed in a particular gauge, it should be understood 
that CTC's conditions themselves are gauge invariant. Also, it may seem quite strange 
that whether the center is well-behaved or not should depend on the starting point of 
the flow. This is probably because the centers in bubble geometries are not singular and 
contain no source. 

We have shown the simplest example when the weak conditions are enough to imply the 
satisfaction of the strong conditions. However, it is unclear if this is true or not in more 
general situations. Nevertheless, the weaker form is apparently far easier to implement 
because we only need to check several points and will be very useful in the numerical 
studies. For the examples we present in section 5, the weak conditions are actually enough. 

With more centers included in the set, the number of possible trees increases very 
quickly^'^. However, the procedure to build the "skeletons" and pick out only those trees 
that satisfy the existence conditions is essentially the same. Therefore, given a set of a 
certain number of centers with some assigned charges and an asymptotic vector, one can 
write a simple computer program to list all of the trees that satisfy at least the weaker 
form of the existence conditions. Furthermore, one can use the "skeletons" of these trees 
to find the corresponding solutions of the bubble equations through assembling the centers 
which will be explained in section 4.3. 

4.3 Disassembling and assembling the solutions 

As we have explained in section 3.3, one can use the skeleton as the specific path for tuning 
the asymptotic vector. By doing so, one can disassembling the corresponding solution. For 
example, one change the asymptotic vector according to the flow depicted in the skeleton 
and continuously traces the deformation of the solution under the constraint of the bubble 
equations. When the asymptotic vector approaches the split point, the solution will divide 



'^'^ Assuming every center is distinguishable, the number of trees is , !• „jv'-i ■ 
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into two independent clusters as the distance between them graduahy increase. After it 
hit the spht point, the system is divided into two subsystems and one continues to fohow 
the spht flows for them individually. Eventually, one can disassemble this solution pairwise 
completely and it serves as an unambiguous way to classify the solution as a particular 
tree. 

On the other hand, the skeletons can also be used as the blueprints to assemble the 
solutions. To find a solution explicitly, we need to fix its sliding degrees of freedom and 
this means to pick a particular point from the whole component. Which point should 
we choose? Although this is not proven, it is quite natural to suspect every split fiow 
component should contain at least one solution that respects the axial symmetry due to 
their nature of pairwise decay. Therefore, for simplicity, we should try to assemble this 
particular axi-symmetric solution. For example, one can try to assemble this particular 
solution for the tree in Figure 3. To do that, we solve the bubble equations of only Fft and 
Tc with /loo set to hg^. Then we bring in the F^ by pushing hoo up flow a little bit and 
probing the axi-symetric solution which has Fq separated very far away from F;, and Tc- 
After we have this solution with three centers, we can change /iqo to ho slowly and trace the 
solution. Of course, for this case with just small number of centers, this procedure seems 
redundant. However, in principle, this procedure can assemble all split flow solutions with 
any number of centers. The only subtleties are the orientations of the different subsets. 
We flnd that sometimes, these orientations are important and one needs to change the 
orientations of some subsets in order to complete the flow. 



5 Numerical examples 

5.1 The merger of two bubbled supertubes 

In this subsection, we review some numerical examples in [9] but look at them from the 
new perspective we have developed. The first numerical example was originally used to 
study the merger of two bubbled supertubes and it turns out to be a good example to 
demonstrate the different aspects of the split flow picture. The basic setup is the following: 
there are two bubbled supertubes and each contains two GH centers with the opposite GH 
charges. Therefore, in total we have five GH centers including the isolated GH center at 
the origin. The charges of a bubbled supertubes can be specified by their GH charge, Q 
and two other parameters: d and X which are related to the charge parameters of the 
individual GH center. If the magnetic part of the charge parameters of the two GH centers 
in the bubbled supertube are: 

F_ = ( — Q ) ^(-) ) ^(-) ) ^(-) ) 1 r+ = (Q ) ^(+) ) ^(+) 1 ^(4-) ) ) (5-1) 

where we assume Q is a positive number and then d^ , X^ are defined as: 

/ = 2(fcf_) + fcf^)) X' ^ 2^-tl-h±l , (5.2) 
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In this example, the charge assignment of the two rings is the following: 



Qa 

Qb 



105, 
105. 



4 



( 50 , 60 , 40 ) , 
( 80 , 50 , 45 ) , 



Xl 



(110,560,50), (5.3) 

(x, 270, 280), (5.4) 



where x was varied from about 64 to the merger value of x ~ 90.3. We are interested in 
axi-symmetric solutions and therefore we put all of five centers on the z-axis and fixed their 
order. After solving the bubble equations numerically, the separations between the center 
can be obtained. Because the number of the center is not so large, we can find all of the 
solutions. In the range of x which goes from about 71.7 to 90.3, we found there are three 
branches of solutions. In the following, we will see how these different branches correspond 
to the different components. We will show two cases with x = 85 and x = 90 in detail. 

5.1.1 x = 85 

For X = 85, three solutions written with the centers' z coordinates are: 





To 


r+a 


r_a 


r+fe 


r-6 


(i) 





1774.222 


1774.225 


2916.575 


2916.581 


(ii) 





2346.036 


2346.037 


2423.872 


2423.873 


(iii) 





2397.797 


2397.816 


2398.288 


2398.318 



The subscript, ita and ±6, indicate they have a positive or negative GH charge and belong 
to the ring a or h. In [9], the solution (iii) has been discovered with a large region of CTC's 
while the other two are free of CTC's. Moreover, the solution (ii) belongs to the branch 
which will eventually reach the merger of the two rings when x ~ 90.3. 

At first, we use the procedure described in section 4.2 to find how many split flow trees 
we have for this setting. For the system with five centers, we have 105 possible trees. After 
going through the weak existence conditions, only two of them still remain. 



rw = {{ro,{r+„,r_4},{r+6,r_,}}, 
r(2) = {ro,{{r+,,r_j,{r+fe,r_fe}}}. 



(5.5) 



To link these two trees to the solutions, at first, we can check the split flow graphs of 
these solutions. We use the normal cylindrical coordinate, {r,9,z) and plot several lines 
defined by (Fj ,'H) = on the (r, z) plane at some fixed 6. The results are shown in Figure 
4. One may notice that they are very similar even for the very different solutions. This 
may be due to the similar underlying structures coming from the same asymptotic vector 
and charge assignments. As the figures contain several irrelevant lines, it is difficult to 
see the pattern of the flows. If we remove the irrelevant and the disconnected lines, we 
have Figure 5. One can see the patterns of the flows become clearer in these graphs. For 
solution (i), the flow starts from F^ and splits to Fo+a and F;, and then Fo+a further splits 
to Fq and F^. Finally, each ring undergoes the splits, F^ (, — )• F_,_(„;,) + F_(-^ ;,), which are 
not visible in the figure. Therefore, we can conclude that solution (i) belongs to the split 
flow component with the tree, T'^'. 
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Figure 4. 

colors are 
(green) , 



From left to right, these graphs correspond to the solutions (i-iii). The lines in different 
drawn by several different conditions: (Tt ,%)={) (blue), (Fq , %) = (red) , (Fa , %) ~ 
(Ff, ,H) =0 (purple) ,(Fo + F^ , "H) =0 (yellow) and (Fa + Ff, , H) =0 (brown). 
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Figure 5. From left to right, these graphs correspond to solutions (i-iii). They are what remains 
after removing the irrelevant lines and disconnected lines. They clearly show the flow patterns of the 
three solutions at least up to the first split. The solution (i) undergoes the first split Tt — ^ Fo+a + F;, 
and then the second split Fo+a -^ To + Fq . The solution (ii) and (Hi) are similar and only their 
first splits, F4 — > Fo + Ta+b, ore visible in the figures. 

For solutions (ii) and (iii), only their first splits, Ft — )• Fq + F^+fo, are visible in the 
graphs. In fact, even we zoom in the near ring region, the next split point cannot be 
recognized due to the bending and disconnection of some lines. To see if the Ta+b flows 
end with scaling solutions or they actually have further split, we can decouple the isolated 
GH center by taking /iqo to the first split point. This means we fix h^o at^^: 



/lo 



0.000417, 0.0271, 0.0229, 0.0177, -3.688, -2.727, -4.152, , 



(5.6) 



and then we solve bubble equations for four centers with this new h. We found there are 
two solutions this time: 





r+a 


r_a 


r+6 


F_fe 


(ii) 





0.00052 


74.8067 


74.8079 


(iii) 





0.0163 


0.416 


0.442 



^^This is the value that we get from the skeleton. One can use the precise value at the split point in the 
split flow graph but the quantitative difference is not significant 
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Figure 6. From left to right, these graphs correspond to solutions (ii-iii). There are several lines 
including: (Ta+'^b ■, "H) = (blue), (Ta , H) ~ (red), {Tb , Ji) = (green) and all other individual 
flows, {T±a,b , H) = (brown). 

Why they are corresponding to solution (ii) and (iii) wih become clear once we study 
the dynamical picture. For now, we can speculate the connection by the ratios of several 
distances. Then the split flow graphs that we get for the solutions are shown in Figure 
6. We can see the first split is clear in solution (ii) but again we have several lines torn 
apart and it is probably because the distance of two rings is too close. However, we can 



proceed similarly and then we see there are indeed further splits: Ta,b — ^ r_| 



ia,b) 



+ T_ 



(a,b)- 



Therefore, we conclude that solution (ii) again belongs to a split flow component but with 
the other tree, T^^'. On the other hand, the graph of solution (iii) has no clear first split. 
This is consistent that it should not belong to a split flow component since the two possible 
trees have been identified. In fact, since there is a large CTC's region surrounding the two 
rings, solution (iii) belongs to a bad component. 

To verify the above identifications, we employ the dynamical picture. As we have 
mentioned, it does not matter which point in the component we trace. For simplicity, we 
focus on the solution which preserves the axial symmetry. In other words, we keep the 
same setup but we change /iqo following the flow described by the skeleton: 



hg^ = mTt + h 



(5.7) 



We start from rji = 0, increase it by a small value at a time and solve the bubble equations 
for each new r/i. The increment should be small enough to trace all of the solutions. The 
first splits of three solutions are shown in Figure 7. There are two split points corresponding 
to the two different split flow trees at: 



vil^ 



(To 



hn 



(ro + ra,r6 



0.00034 . 



vif 



( To , /io 



( Fo , Fq + Fft ) 



0.00042. (5.8) 



The subscript "si" denotes the first split and the superscript label the tree it belongs to. 
We can see these two split points clearly in the figure. When the value of rj is close to rjsi , 
the two rings in solution (i) indicated by the top and the bottom blue lines are quickly 
separated far apart and after rji crosses r/^j , solution (i) disappear and we left with only 
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Figure 7. This graph traces the z-coordinate of all centers for all solutions under the change of rji. 
As there is not much change when 771 = '^ 0.0002, we show the graph starting from rji = 0.0002. 
The red line along z — represents Tq for all solutions. The positions of the other centers are 
shown with blue dots. The remaining centers of solutions (ii) and (Hi) are very close to each others 
and are represented by the middle blue line. The top blue line is Ft, and the bottom blue line is 
Ta for solution (i). The brown dashed line marks the exact value of rfs^ and the green dashed line 
marks rjs^ . 



two solutions. Notice that all of the rings in solutions (ii) and (iii) are very close to each 
other and accumulate at the middle blue line. When ?/i increases further and is close to 

(2) 

Tjsi , the separation between rings and Fq for the remaining two solutions begins to increase 

(2) 

very rapidly. After r]i passes r/^/, no solution exists. This analysis confirms the existence 
of the first split for all of the solutions. 

For the second split, we focus on solutions (ii) and (iii) because the split-flow pattern 
of solution (i) is already quite clear in Figure 5. We consider that hc^ starts from r]i = j/^i 
and follows the flow of Ta+b'- 



(2) 



hg^ = r/2(F„ + F6) + hs, = r?2(Fa + Ffe) + lif Ft + h. 



(5.9) 



where we parametrize this flow with 772 and start the flow from 772 = 0. Then we solve the 
bubble equations for four centers and the result is shown in Figure 8. The next split point 
of T*-^^ can be easily computed: 



„(2) 
IS2 



(r, 



a ) "-si / 



(r« 



0.0167. 



(5.10) 



We can clearly see that for solution (ii), Fq and F;, are separated further away when r/ 

(2) (2) 

approach r/^j and after 772 crosses r/^j , the solution disappears. The following split points 

of solution (ii) can also be clearly identifled and we again conclude that it belongs to T^"^' ■ 

On the other hand, the behavior of solution (iii) is completely different. The overall size of 

(2) 
solution (iii) decreases as T72 increases and even after 7^2 crosses 77^2 , it continues to exist 

with decreasing size. This is the typical behavior of the scaling solution and this part of 
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0.02 0.04 



Figure 8. From left to right, these graphs correspond to solutions (ii-iii) and they show how the 
solutions change under the flow of 7)2 ■ The top (double) blue indicates Ff, while the bottom (double) 

(2) 

blue line represents Ta- The red dashed line marks the exact value of rjs^ . 

solution (iii) should be regarded as a single flow from the point view of the split flow picture 
. However, as we have already mentioned, solution (iii) has CTC's and should belong to a 
bad component. This concludes that the dynamic picture is consistent with the static one 
in this case. 

5.1.2 x = 90 

For X = 90, we still have three solutions: 





To 


r+a 


r_a 


r+6 


r_6 


(i) 





1774.10378 


1774.10683 


3005.96801 


3005.97352 


(ii) 





2419.88471 


2419.88474 


2423.38597 


2423.38604 


(iii) 





2422.25748 


2422.25842 


2422.27682 


2422.27830 



However, if we again look at the physical trees, we find that only one remains: 7"'^^. The 
reason that T^"^' vanishes is because (r^ , Tf,) change sign and consequently r/^j become 
negative. To understand why we still have three solutions, we again employ the static and 
dynamic analysis. We find that the behavior of solution (i) does not change much and it 
still match with T^^' ■ On the other hand, solution (iii) still has CTC's and should belong 
to a bad component. Therefore, we focus on solution (ii). Furthermore, since the first split 
point still exists just like the previous case, we focus on the analysis after the isolated GH 
center is decoupled. Proceeding similarly, we have the result shown in Figure 9. We see a 
curious fact that even though the solution itself does not change much, its behavior under 
the flow dramatically changes. From Figure 9, it is clear that now solution (ii) shows the 
scaling behavior. Since solution (ii) is free of CTC's unlike solution (iii), it truly belongs 
to a scaling component. To put it more precisely, the T^^^ vanishes at exactly x = ^^ 
when {Ta , F^) flips the sign and that is exactly^^ the point when solution (ii) begins to 
show the scaling behavior. This concludes that when ^P < 2; < 90.3, solution (ii) belongs 
to a scaling component. 



^^For the case with x only shghtly smaller than the critical value, ^^1^, it will be difficult to determine 
the split point in the split flow graph. However, the solution continues to show the splitting behavior in 
the dynamic analysis. 
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Figure 9. On the left, we have the split flow graph for solution (ii) after the first split. The blue 
line is (Fq + Ff,, H) ~ and the green and red lines are {Ta^b, "H) = 0. On the right, the graph 
shows how solution (ii) changes under the flow with the top blue line indicating the positions of Tb 
and the bottom blue line marking the positions of Ta ■ 

Another interesting configuration is x ~ 71.7 when solutions (i) and (ii) become very 
similar. However, there are still two trees when —^ > x> 14.156^^. If these two solutions 
somehow merge to one solution, the split flow picture will not be correct. However, inter- 
estingly, no matter how much we fine tune the value of x, we always get either both of the 
solutions or neither of them. With the help of flux quantization, these two components 
may be separated without touching each other. Also, it may seems peculiar why both of 
them vanish when x < 71.7 while the corresponding trees still exist. The reason is simply 
that the order of the centers sometimes becomes important and when x < 71.7 the solu- 
tions corresponding to these two trees exist with the different order of the centers or more 
general orientations without axial symmetry. 

5.2 Triangular solution 

In [9], it was found that if the fluxes between the three centers form a closed triangle, 
the scaling solution exists. Specifically, one can adjust the angles between three centers to 
make the distances between them shrink and eventually merge into a point in the coordinate 
space. Physically, this produces a very long throat and makes the solution become very 
similar to a classical black ring with a finite horizon area when the three centers are very 
close to each others. Also, an example of this kind of scaling solution was given. We would 
like to see if this example shows the scaling behavior from the perspective of the split flows. 
The setup is the following: 







(1,0,0,0), 

(-105, 1200, 16000, 1613), 



^ 1 

pM 

^3 



(-105, 2210, 7887, 7800) , 
(210, 525, -20000, 6400) . 



For this setting, there is an axisymmetric solution. However, no tree exists for this set 
^^When X < 14.156, (Ta, hs-^) changes sign and again the risj becomes negative and T'^' vanishes. 
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Figure 10. On the left, the graph show how the triangular solution changes under the flow before 
the split. The red line indicates Fq while the other centers accumulate at the blue line. The yellow 
dashed line marks the critical value, rjs ~ 0.00172, for the split. On the right, the graph traces the 
solution with the remaining three centers afler the split. 

of charges. This imphes that this solution should belong to either a scaling or a bad 
component. Since it has been verified that no pathology is present, this solution indeed 
should belong to a scaling component. We can implement the dynamic picture to confirm 



(m) T^{m) 



^(m) 



this. Since the scaling solution is formed only by F^ : ^2 ^^^ ^3 ' ^^^ split between 
them and Fq still exists. This split is shown clearly in Figure 10. After the isolated GH 
center is decoupled, we continue the fiow for the three centers in Figure 10. The scaling 
behavior appears at this stage and one finds that it cannot be torn apart by tuning the 
asymptotic vector along the flow. This confirms that this triangular solution indeed belongs 
to a scaling component. 

5.3 Three rings 

To increase our confidence in this picture, we would like to study more complicated solu- 
tions. For example, we can have three rings and this means in total we have seven centers. 
We take the charge assignment of the three rings as the following: 



:50,70,65), Xi 



Qa = 105, di = 

Qb = 105 , dl = { 106 , 55 , 35 ) , XJ^ 

Qc = 105 , 4 = ( 35 , 80 , 100 ) , X^ 



;920, 880, -270) 
; 310 , 590 , 600 ) , 
; 520 , 880 , 220 ) . 



(5.11) 



For this case, we find that only three trees pass the weaker existence conditions: 

r« ^ {{Fo,{r+,,r_4},{{F+fe,F„fe},{r+c,r_j}}, 
T('^ ^ {{Fo,{{F+„,F_j,{r+,,r_j}},{r+,,r_a}, 
T^'^ = {{{ro,{r+„,F_j},{F+,,F_j},{r+b,r„f,}}. 

For this many centers, one will have difficulty finding all of the solutions even numerically. 
However, by the method described in section 4.3, one can assemble all of the solutions 
corresponding to these trees. For simplicity, we still look for the axisymmetric solutions. 
The corresponding solutions we find are the following: 
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Figure 11. From left to right are the split flow graphs for solutions (1) to (3). The red line is 
(Fq, %) = 0, the green line is (Ta, ^) = 0, the yellow line is (Vh, T~L) = and the purple line is 
(Fc, "H) = 0. The other color lines are for some combinations of the charges. The splits between 
two centers in the rings are not visible in the graphs. 





To 


r_a 


r+a 


r-6 


r+fe 


r-c 


r+c 


(1) 





7166.711 


7166.712 


16693.003 


16693.006 


20337.983 


20337.987 


(2) 





13615.977 


13615.979 


20481.737 


20481.739 


11916.918 


11916.922 


(3) 





8654.738 


8654.741 


-18007.984 


-18007.988 


16105.593 


16105.598 



As usual, we express the solution with the z-coordinate of each center and fix Fq at the 
origin. Notice that for solution (2) and (3), we need to rearrange the order of the centers 
in order to find the solutions. Their split flow graphs are shown in Figure 11. 
Additionally, we find another axisymmetric solution for this set of charges: 





To 


r_a 


r+a 


r-i 


r+6 


r_c 


r+c 


(s) 





14228.9671 


14228.9676 


14809.0745 


14809.0752 


14557.1942 


14557.1951 



This solution does not belong to any tree and in fact, should be classified as a scaling 
solution. To see this clearly, one can bring the asymptotic vector to the only split point 
and decouple Fq. One can then carry out the dynamic analysis for the remaining six 
centers. The result is shown in Figure 12. Their scaling behavior is quite clear. 

The structure of this scaling solution is actually similar to the triangle solution intro- 
duced in previous subsection. One can see that the fluxes between three rings indeed form 
a closed triangle. Therefore, the existence of this solution implies that the scaling solutions 
may have multi-layer structures just like the split flow solutions. It also show that the 
scaling and split flow components can coexist and a set of centers can have multiple split 
flow and scaling components at the same time. 

5.4 Bifurcation points: splittings of the components 

In this subsection, we will show the possibility that a connected component can be split 
into two components under the flow. We will use the example in section 5.3 with only some 
minor changes. We change only d\. to the following and keep the other settings the same 
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Figure 12. This graph traces three rings under the flow after Tq is decoupled. The green line is 
Ta, the yellow line is Ff, and the purple line is Te- 



as (5.11): 



di = (35, 20, 100) . 



(5.12) 



For this setting, we have the similar scahng solution but different trees that pass the weak 
existence conditions: 

r(^) = {ro,{{r+,,r_,},{{r+„,r_4,{r+e,r_j}}}, 
T^'^ = {{ro,{{r+„,r_j,{r+,,r_j}},{r+b,r_b}}, 
T^'^ ^ {{{ro,{r+„,r_j},{r+fe,r_5}},{r+c,r„j}. 

Among these three, T^^' is the most interesting because it has the same first split point as 
the scaling solution. This means that after we decouple Fq, we should have two components 
corresponding to T^^' and the scaling solution respectively. However, we find that there is 
only one component when we take /iqo to the split point: 



h 



Sl 



^tVsi + ho, 



(5.13) 



where Ff is the total charge and r/^j is equal to: 

(ro, ho) 



Vs^ 



(ro, Tt) 



0.000064805 . 



(5.14) 



If we change hoo to continue the flow: 



ha 



{Ta + Tb + T,)r]2 + hs. 



(5.15) 



we find that at some point, the component breaks into two parts. We show the axisym- 
metric solutions in the following table: 
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T-a 


r+a 


r_b 


r+6 


r_e 


r+c 


??2 = 





0.0021 


1454.3822 


1454.3830 


8066.9990 


8067.0021 


V2=r]*2, (i) 





0.00158 


12269.368 


12269.370 


-1280.6451 


-1280.6444 


m = vh (ii) 





0.00158 


12271.433 


12271.435 


-1280.6819 


-1280.6812 


7?2 ~ 2.65 X 10-^ (i) 





0.00147 


7186.637 


7186.639 


-1134.0126 


-1134.0119 


rj2 ~ 2.65 X 10"^, (ii) 





0.00165 


29530.173 


29530.176 


-1409.4583 


-1409.4576 



From the table, one can see that at the beginning of the flow (772 = 0), only one component 
exists. However, starting from ?72 ~ ^2 ~ 2.5459751757 x 10~^, we have two components, 
solution (i) and (ii). 

To understand what is happening, we need to study the whole solution space. In this 
situation, this space can be effectively captured by a three-center system. In other words, 
we will analyze the solution space of the coarse grained system by replacing each ring with 
only one center. For 772 = 0, the independent part of the bubble equations can be written 
as: 



Tab 



rgb 

nc 



{Tb,hs,). 



(5.16) 
(5.17) 



As {Tg, Tc), (Ffe, Tg) and (Fc, F^) are all positive and satisfy triangular inequality, the scal- 
ing solution of these three centers exists. Furthermore, there are at most two components 
in the solution space (for more detail, see appendix A in [21]). To know how the compo- 
nents change along the flow, we can parametrize the component with one parameter. A, as 
the following: 



rgb 



A 



rbc 



{rb,ra) 



(rcTft) ^ y 



A 



(Fa.rc) 



+ 



(5.18) 



where x, y and z are the three distances of a particular solution in the component. For 
example, one can use axisymmteric one as this particular solution. Then, for any value 
of A such that all distances are positive and the triangular inequalities are satisfied will 
correspond to a solution. Therefore, the solution space in this case is the collection of 
some (at most two) intervals on a one dimensional line parametrized by A. To check how 
this space changes along the flow, we take the axisymmetric solution with y + z = x and 
monitor rgj, + rcg — rbc for several points along the flow. The result is shown in Figure 13. 
For A < 0, either some distances will become negative or some triangular inequality (i.e., 
rbc + rcg > rgb) will be violated and the only possible domain is positive A. In Figure 13, 
the valid region must have the positive value of rgb + rcg — rbc, it is clear that the interval 
(i.e., [0,00)) splits into two intervals (i.e., [0,A_] and [A+,00)) when r/2 is larger than 772- 
From (5.18), it is clear that A — t- 00 corresponds to the scaling point. Therefore, we clearly 
see that along the flow, the scaling component splits into two in which one of them becomes 
the split flow component, T^^' ■ 
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Figure 13. This figure shows how rat + rca — ri,c changes with A for several different values of rj2- 
The blue curve corresponds to 772 = 0, the red curve corresponds to 772 = rj2 and the green curve 
corresponds to 772 = 2.65 x 10~^. When one increases the value 0/772, the curve becomes gradually 
lower and eventually, some part of it becomes negative which violates the triangular inequality. 

Since we can continuously change the background, /iqo, there should be a precise bi- 
furcation point when the component just splits into two. This precise value of 172 for 
this bifurcation point naturally must be smaller than the next split point of T"'^-* (i-e., 
r]s2 ~ 2.9 X 10~^). From the point view of the split attractor flow, when the flow passes the 
wall of marginal stability, it can split or continue the flow toward the single attractor point 
if this point is well-defined. What we show here indicates the system must decide its fate 
precisely at this bifurcation point when the component just splits into two. Interestingly, 
at this moment, if the solution is exactly at the joint point of the two components, it 
cannot decide whether it should fall into the scaling component and have all of the centers 
merging into the scaling point or it should fall into the split flow component and completely 
reduce to several decoupled GH centers. The physical implication of this is unclear at this 
moment. However, one thing is clear that one cannot distinguish the split flow from the 
scaling components unambiguously before the appearance of this bifurcation point. 

6 The existence of three-center trees 

In bubbled geometries, the solutions are composed of several centers and each center is 
basically parametrized by a three dimensional vector, fe[, and its GH charge. In order to 
have solutions which are free of CTC's, all of these charge parameters must be constrained 
in some way. How they should be constrained in general is not clear. However, for the 
solutions that belongs to split flow components, the existence conditions of the tree should 
immediately provide these constraints assuming the conjecture is true. In this section, we 
will see how these constraints emerge from the existence conditions by using the next-to- 
simplest tree with three centers as an example. 

For simplicity, we assume every center has a GH charge of either -|-1 or —1. Also, 
we focus on the case in which the metric is asymptotically flat and this means that the 
GH charges sum to one and hoc is always set to /iq defined in (2.22). Given a set of 
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three centers, (Fq, Ta, r^), the possible trees are {Tq, {Ta, ^b}} with two other cychc ones 
(see Figure 14). Naively, we will have nine parameters coming from the remaining charge 
parameters but we can remove three of them by using the gauge symmetry in (2.26). We 
can always choose to fix the gauge such that the isolated charge becomes trivial and the 
resulting charges are: 

F^'") = ±(1,0,0,0), ri-) = (1, kl kl kl), Ff ) =^(1, kl kl kl), (6.1) 

If the isolated charge is positive, the solution will have the same split flow pattern as the 
simplest bubbled supertubes. Therefore, it is convenient to implement the ring charge 
parameterization and define F^ as: 

F, = Ta + n = (0, p\/,/^-sipV, -S2pV, -S3pV, 2s„pV/), (6.2) 
where {p^-,sj) are related to {ka,kh) by: 

^a ~ H — P 1 kg^k^ — kfj kjj = — sip p , 
klkl-klkl = -S2p'p\ klkl-klkl = -S3P^p\ (6.3) 

and Sm is not free but fixed by sj: 

Sm = -{I- sj- S^- sl + 2s2S3 + 2siS3+2siS2). (6.4) 

o 

Notice that sj and Sm are all dimensionless quantities. To further discuss the constraints, 
we will divide the analysis into three sectors: 

• Sector (I): The isolated charge is positive (i.e., Fq = (1,0,0,0)) and p^ are either 
all positive or all negative. 

• Sector (II): The isolated charge is positive and one or two of p are negative while 
the remaining are positve. 

• Sector (III): The isolated charge is negative (i.e., Fq = —(1,0,0,0)). 

These three sectors cover all possibilities for the three-center trees in bubbled geometries. 
We will analyze them in the following subsections. 

6.1 Sector (I) 

In this sector, Fq is positive and p^ are either all positive or all negative. The tree is 
T = {Fq , {Fa , Ffo}} and there are in total five edges in the tree as shown in Figure 14. 
The strong existence conditions of the tree require there is no CTC along any of these 
edges. However, we will show that the weak conditions are enough. In fact, the sufficient 
conditions are the even smaller subset of that for this sector: 

• The separation of the first split is positive: Vsi = :f- > 0. 

• The separation of the second split is positive: r/^j = r — ^ ^' 
. VZi{ro)>0, 1 = 1,2,3. 

The first condition is: 

%i = — = - ir T' \ = n r^-3 > 0- 6.5) 

rsi (Fo, Fr) 2smP^P^P'^ 
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(I) 




Figure 14. The structure of a general 3-center tree. There are in total five edges labeled from 
(TV). There are two split points labeled with ha^ and hg^. 



In this sector, this imphes Sm must be positive. The second condition is: 

1 _ (Fa, /isi) _ (r^, To +rfe) + (Fa, /loo)rsi 



^s 



where hg-^ is: 



(Fa, Tfe) 



h. 



plp2p3 



Vs^ 



[To + r,.) ?7si + he 



and A^ is defined as: 

As = 2Sm(l -/liSi -/i2S2 - A^SSs) - l + Sl'S2S3, 



AsTjs-,, (6.6) 
(6.7) 

(6.8) 



with fij = ^ — J. The positivity of rjs^ and r/^^ implies that A^ must be positive. The 
last condition requires VZi to be positive at Tq. As we have explained in section 4.2, this 
is equivalent to the conditions in (4.8). Using the explicit form of /i^^ in (6.7), we have the 
following conditions: 



1 + ^i^ > 
1 + ^^^ > 



1 2 



'SI 



> 



2/UiSm + Si > 0, 
2/U2Sm + S2 > 0, 
2H3Sm + S3 > 0. 



(6.9) 

(6.10) 
(6.11) 
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To summarize, we have the following conditions: s^ > , A^ > and (6.9-6.11) in this 
sector. These conditions impose the constraints on the charges. However, before we can 
claim these constraints define the physical region, we need to verify that the strong existence 
conditions of the trees are satisfied. Due to the complexity of the CTC's conditions, we 
start with the parallel charge case to get a better idea about the basic features of these 
constraints and then move on to the general case. Additionally, we also check the dangerous 
half-split point to verify the conjecture. 

6.1.1 Parallel charges 

For the parallel charge case, we simply set s/ to be equal (i.e., si = S2 = s^ = s). 
Immediately, we have Sm = g(l + 3s^) and the positivity of it is given. The remaining 
conditions become: 

s^ + 3s^-s-3 > 0, 3s'^ + 12s + l > 0, (6.12) 

where the second inequality comes from the summation of the conditions (6.9-6.11). Com- 
bining the above two conditions yields s > 1. Therefore, s > 1 defines the physical region 
for the parallel charge case. This also implies that the total charges of harmonic function 
Li must be positive. Consequently, Lj is positive along edges (I-III) and VZ^ > is triv- 
ially satisfied along these edges. In the following, we will check the remaining unverified 
CTC's conditions for all of the five edges in the physical region. 

Edge (I): Since VZj > is always true for s > 1, we only need to check the other 
condition, Q > 0, along this edge. To monitor how Q changes along the flow, we expand 
Q with the power of ^^: 

Vsi 

and the coefficients are: 



or 


= , qi^^ = r]si , qf ^ 


= {p^p^ + p^p 


qf' 


I 




qS" 


= - ^ (pW)' (1 + sf (17 - 


- Ws + 9^2) ,?4^ 



(6.14) 

For s > 1, q^ is always negative while the other coefficients are all positive. This indicates 
that if we put all of the charges into a single center, it will not be a well-defined attractor 
point and implies the single center solution or the scaling solution does not exist and the 
first split must happen. To check Q remains positive along edge (I) even with the negative 
q^ , we will show that q^ > 4|qj^ | in the physical region. For the parallel charge case, 
this is immediately followed by: 

(l + Ss^) > -(17-10s + 9s2). (6.15) 
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This inequality is always true for s > 1 and saturates at the boundary, s = 1. Because of 
qi > 4|q]^ I, only the third and fourth order terms are enough to guarantee the positivity 
of Q at the split point where ri\ = r/^^. In fact, this is still true even when we go deeper 
into the half-split point where 771 = 2?7s^. Curiously, we see that the weak conditions are 
enough to ensure Q > at this dangerous point as well and avoid the conjecture to be 
violated. Furthermore, if we take the derivative of Q with the dimensionless parameter, 
^, we see that the derivative is always positive along edge (I) (actually, it is still positive 
up to the half-split point). This means that Q monotonically increases with r/i for this 
whole edge. 



Edge (II): We can do a similar expansion of Q with ^^ for this edge. The coefficients 



"&' = pVpM/+L^)3 ((^')^^^^(-) + {V^')i''\s)),.,, (6.16) 



(2) = /"f-pSwCS). X , (vIJ\d^)(,\\„'^ 



q^ = [{V')e'Ks) + {V'')^''^{s)Ul, (6.17) 



(3) P VV El P (g +3s'-S-3) 3 , 1 2 3n2 4 /«in^ 

q2 = ^^3^2) ^«2 = (PPP) Vs2^ (6-18) 



(4) _ (plp2p3^2 

4 '"^2 



qr = - ^" Y ' Vi , (6.19) 

where V"^ and V^'^ are defined as 



V"^ = p^p^ +p^p^ +p^p^ , V^-^ = p^p^ +p^p^ +p^p^ , (6.20) 

and ^(s) are some polynomials of s: 



c«(.) = 


- 5 + 385^ + 485^ + 15s^ , 


(6.21) 


c^'Hs) = 


- 8 + s + 885^ + 102s^ + 48s^ + 9 s^ , 


(6.22) 


C^'\s) = 


- 13 - 6s + 34s2 + 305^ + 3s^ , 


(6.23) 


&Hs) = 


- 25 - 12s + 7is'^ + 60s^ + 15/ . 


(6.24) 



One can easily show that these four polynomials are all positive when s > 1. Therefore, 
q2 and qg are positive. Furthermore, one notices that in this edge, q^ is exactly equal 
to 4 |q2 |. This tells us that not only is Q positive but also that it increases monotonically 
along this edge. Additionally, it also shows that the dangerous half-split point (i.e. 772 = 
27/^2) for this edge is also safe with positive Q. 

Edges (III-V): All of these edges are single-center fiows that terminate at some GH 
centers. As we have mentioned, we only need to check the conditions at the start and end 
points for these kinds of flows. All of the starting points should have been verified in the 
previous two edges and therefore, we only need to check the conditions (4.8) to ensure the 
end points are also well-behaved. 
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For edge (III), we have already required the conditions in (6.9-6.11). For the remaining 
two edges, we start with hg^: 

hs, = ^J^±^ + hs, . (6.25) 

To apply the condition to Fq and F;,, we change the gauge respectively such that Fq or F;, 
explicitly become pure a positive/negative unit GH center and then check the electric part 
of hg^ . The conditions for F^ are: 

sil + s"^) l + 5s2 ^ ^ 

+ , : , o > /= 1,2,3, (6.26) 



and for F;, are: 



/i/(l + 3s2) l + 3s2 



■ 1 > 1 = 1,2,3, (6.27) 

///(l + 3s2) 

where ^/ = yfi — f and since p^ are all positive or all negative, fij are between and 1. It 
is clear that these conditions are satisfied for s > 1. Therefore, for the parallel charges, we 
conclude that s > 1 is the only necessary and sufficient condition in order for the trees to 
exist and it defines the physical region. 

6.1.2 General charges 

In this subsection, we remove the restriction on sj. Just like the parallel charge case, we 
take p as some arbitrary constants which are either all positive or all negative and check 
the constraints imposed on s/ by the sufficient conditions. At first, we divide the three 
dimensional s/-space into eight regions according to the signs of sj (i.e., (ib,ib,zb)) and 
look at the conditions for each region. Interestingly, we find that no matter what the values 
of p^ are, the physical region is always confined in the all-positive region (i.e., (+,+,+)). 
In the following, we will assume si > S2 > S3 but a similar argument can apply to the 
other cases by taking advantage of the cyclic nature of the conditions. To start, we express 
As as the following: 

As = As + 2 ((S3 -Si)^i + (S3 -S2)/i2) Sm, (6.28) 

where As is independent of fii and is defined as: 

As = 2Sm{l-S3) - (I-S1S2S3), (6.29) 

also we have used /i3 = 1 — ^i — /.i2 and Sm is defined in (6.4). At first, we check the regions 
except the all-postive region: 

• (si, 82,53) G (+,+,—): Because Sm must be positive, the term that depends on fij 
in (6.28) is always negative. Therefore the maximum of As is simply As. However, 
one can show that As is negative in this region. Therefore, it is clear that As and 
Sm cannot be simultaneously positive in this region as long as p^ are all positive or 
all negative. 
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• {si, 82,83) G (+,—,—): To exclude this region, we need to include the other three 
conditions (6.9-6.11). By combing (6.10) and (6.11), we have: 

2Sm + S2 + S3 > 2lIiSm- (6.30) 

Because fii and Sm are positive values, this implies: 

2Sm + S2 + 83 > 0. (6.31) 

On the other hand, to maximize A^, we can set /xi to zero but not fi2- Because S2 
is now negative, we can only take the minimum value of fi2 allowed by the condition 
(6.10). Therefore the maximum value of A^ is: 

As{fii = 0,fi2 = -^) = A, + 82(82-83). (6.32) 

One can show this maximum value of A^ is always negative under the condition (6.31) 
in this region. Therefore, this region is also excluded. 

• (si, 82, 83) G (— , — , — ): At first, by combing all of the three conditions (6.9-6.11), we 
have the following: 

Sl + 82 + 83+2Sm > 0. (6.33) 

To maximize A^, we need to take the minimum values of /zi and //2 allowed by the 
conditions (6.9) and (6.10) and we have: 

As{fii = -- ,fi2 = -i: ) = As + 82 {s2 - 83) + si {si - S3) . (6.34) 

Z Sm Z Sm 

One can show this maximum value of A^ is always negative under the condition (6.33) 
in this region. Therefore, this region is also excluded. 

Through this analysis, one can conclude that the physical region must be confined in the 
all-positive region. Since the conditions (6.9-6.11) are trivially satisfied in this region, the 
precise physical region is defined only by the positivity of A^, Sm and s/. Of course, 
the exact physical region in s/ space will depend on the value of p^ or more precisely ///. 
However, the general feature of this region is roughly the same with the boundary changing 
a little with /x/. For the special case with the equal dipole charges, the physical region is 
shown in Figure 15. 

It is quite interesting how these conditions interplay and the constraint that s/ must 
be all positive emerges. Similarly, this implies the harmonic functions, L/, are all positive 
along edges (I-III) and VZj are also positive along these edges. In the following, we will 
check the remaining CTC's conditions along all edges in the physical region defined by the 
positivity of Sm, As and s/. 
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Figure 15. The physical region defined by As > and s„i > in the positive sj region with 
Ml = m = M3 = I • 

Edge (I): The coefficients of the expansion of Q with ^ in the general case are: 



(1) 

qi = Vsi 

.(3) 



qi 



(2) 



(P2P3(l + 'Sl) +PlP2(l + S2) +PlP2(l + S3)) ril^ , 



q^3J = ^y^s (pl(l + s2)(l + S3)+p2(l + Sl)(l + S3)+/(l + Sl)(l + S2)) r,l,, 
q(4) = (/pV)2 ('(l + s,)(l + s2)(l + s3)-i(2 + 2s^ + si+S2 + S3)2') <(.6.35) 

For the positive s/, only q^ can be negative. However, we will again show that q^ +4qi 
is always positive in the physical region. 

In the following, we will assume S3 > S2 > si but a similar argument can apply to the 
other cases by using the cyclic nature of the conditions. At first, we isolate the /i dependent 

qf)+4q;^) = (pW)2<qI, (6.36) 

(6.37) 



terms from q^ + 4 q^ as: 



with qi is defined as: 

qT = q! + Ml(l + S2)(S3 -si) + M2(l + si)(s3 - S2), 
where qj: 

ql = 4(l + si)(l + S2)(l + S3)- (2 + 2 s„ + si + S2 + 53)^ + 2 s„(l + si)(l + S2), (6.38) 
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where //3 = 1 — /ii — /i2 is used. To verify q\ + 4qj^ > 0, we will show that the minimum 
of qi is positive in the physical region. As we assume S3 is larger than the other two 
parameters, we need to set //i and //2 to their smallest allowed values to minimize qi". On 
the other hand, the parameters, fii and /X2, are constrained by the positivity of A^: 

N 1 ~ Si S2 S3 

fJ-i [S3 - si) + 1J.2 (S3 - S2) > 1 + S3 . (6.39) 

Z Sin 

The minimum of qT is either at /ii = or ^2 = or /xi = /X2 = depending on the value 

of s/. For the case when both of them are zero, one can show that qi is positive in the 

physical region. For the case with fii = 0, we can take ^2 to its minimum allowed value by 

(6.39): 

I-S1S2S3 I-S3 

^^2m^n = 7, 7 7 ' <^"-'^0) 

^Sm[S3-S2) S3 - S2 

and in order for this to be a meaningful minimum, it must be between zero and one. Under 
this constraint, one can show that the minimum of qi is still positive in the physical region 
for this value of ^2- Similarly, this is also true for the case where ^2 = 0. Therefore, this 
concludes that qj is always positive in the physical region. As we have explained in section 
6.1.1, this also implies that Q increases monotonically with rji along the whole edge (I). 

Edge (II): The coefficients of the expansion of Q with ^ in this edge are: 

^SmPlP2P3 ^4s^ 4Sm 

+ Ail M2 (S3 + 2) + fl2 fJ'3 {si +2)+ fli fl3 (S2 + 2) + 1 ) r]s^ , 



(2) {P1+P2+P3) /^ 3 , 

q2 = 5 (2A3 + I) +^(s/,^7) + 8Sm(//l/X2 + /i2/^3 + /il^3) W/, 



8s 



m 



ls2 ' 



qf = {p'pVf^li, q?^ = -\{pVp'?r^i, (6-41) 



where ^(s/,/x/) is: 



Cisi,IJ,l) = 2(si + S2 + S3) +4(s2S3/ii + SiS3^2 + SlS2^3) -SSm. 



(6.42) 



Notice that the positivity of qg should already be verified at edge (I) and we do not need 
to check it here. From (6.41), we immediately see that q2 is always negative. However, it 
is also clear that qg + 4q2 = 0. This means that q^ is large enough to overcome q^ 
before this flow reaches the split point at 772 = ??S2 01^ even deeper to the half-split point 
at rj2 = 2rjs2- For the positivity of the remaining two coefficients, we only need to check 
^(s/,/i/). With similar arguments as we have been using, one can show ^(s/,/z/) is always 
positive in the physical region. Therefore, we conclude that not only is Q positive along 
this edge but it also increases monotonically with rj2- 

Edges (III-V): For these three edges, we need to check the conditions for single 
center flow. The condition for edge (III) has been included in the definition of the physical 
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sector and is automatically satisfied. For the remaining two edges, we just need to generalize 
the conditions in (6.26, 6.27). The conditions for Fq are: 

1 + ^^ (4A, + (s2-S3)2-s? + 2si+3) > 0, (6.43) 

1 + -^ (4A, + (si-S3)2-s^ + 2s2 + 3) > 0, (6.44) 

1 + ^^ (4A, + (si-S2)2-si + 2s3 + 3) > 0, (6.45) 



8fl3Sr. 



and for Th'- 



-1 + ^^ (4A, -(s2-S3)2 + s? + 2si + 5) > 0, (6.46) 

-1 + -^ (4A, -(si-S3)2 + si + 2s2 + 5) > 0, (6.47) 

-1 + T^ (4A, -(si-S2)^ + si + 2s3 + 5) > 0. (6.48) 

SflsSm 

Using similar arguments, one can show that all of these conditions are satisfied in the 
physical region. 

This concludes the verification that the strong existence conditions of the trees are 
satisfied in the physical region for this sector. According to the conjecture, every existing 
tree should correspond to some well-behaved supergravity solution. In other words, there 
exists a well-behaved solution for every point in the physical region (e.g. Figure 15) with 
arbitrary all-positive or all-negative dipole charges. Before we can claim these solutions 
are microstate geometries of some zero-entropy black rings, we need to verify the following 

ratio: 

(F, + Fb , Fo) 



,r r\ 2s^. (6.49) 

This ratio measures the relative size of the bubble that supports the ring's radius and the 
bubble that nucleated inside the ring. For a three-center solution in bubbled geometries to 
resemble a classical zero entropy black ring, this ratio should be much larger than one (i.e., 
z^ ^ 0{1)). Since this ratio is just proportional to Sm, one can see that it is not bounded 
above in the physical region (Figure 15). This can also be easily seen in the parallel charge 
case. The only dangerous part is near the boundary where Sm goes to zero. When the 
solution is close to such a boundary, the ratio can become comparable or even small than 
order one and it will no longer resemble a black ring. However, the majority of the solutions 
in the physical region still have a large ratio and can be regarded as the simplest bubbled 
supertubes. 

6.2 Sector (II) 

For this sector, one or two of the dipole charges of the ring are negative. Naively, it seems 
unlikely to have well-behaved solutions in this sector. However, we will show that the trees 
still exist in some limited region. 
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To simplify the analysis, we reduce one degree of freedom and take si = S2 = s. We 
will use the following parameterization: 

r(™) = {0,p\p^p^) = {0,p + q,p-q,-yp), (6.50) 

where y is a dimensionless quantity. For this sector, we will focus on the region where p^ 
and p^ are positive while p^ is negative and without loss of generality, we assume p^ > p^ 
(i.e., p > g' > and y > 0). Notice that all of CTC's conditions are invariant if we flip the 
signs of all k^ simultaneously (of course, the asymptotic vector must change accordingly). 
Therefore, the case where two p are negative should be the mirror of the case we analyze 
here. 

We will follow a similar procedure to the previous section. Due to the simplification, 
the sj--space now is only two dimensional and we divide it into four regions based on the 
signs of s and S3. By carefully checking the weak existence condition (particularly, CTC's 
conditions for Fq, F^ edges and the positivity of the two separations.), we find that the 
physical region is confined in (s, S3) G (—,+). For this region, the charges of the harmonic 
functions, Lj, are again all positive. Furthermore, the only conditions we need for this 
region are the ones that require two separations are positive: 

4(j/-2) ^ p ^ (^ 

(1 + 4SS3-S2)(p2_g2)y (l + 4ss3_52) 



^^1 = ,.,.„„ 2N ; 2 2^ > =^ 7V^ —^, > 0' (6-51) 



{l + Ass3-sl){p'^-q^)y (l + 4ss3-s^) 



^^^ = 7^^-T .2w„2_.2u, > =^ n^...„_.2, > 0' (6-52) 



where ^ is: 
ay,s,S3) = (l + 4ss3-si)f(?/-2)(l-S3) + 2(s-S3)) + 4(s2s3-l)(y-2). (6.53) 



Notice that when y = 2, the three dipole charges sum to zero and the separation of F,. 
and Fq diverge. This means that the /iq is already at the split point and consequently, 
the tree cannot exist. For < y < 2 and y > 2, the physical region where the tree exists 
is bounded by two lines in (s, S3)-space. These two lines are defined by the following 
conditions respectively: 

l + 4ss3-si = 0, ^{y,s,S3) = 0. (6.54) 

For < y < 2, the first line provides the lower bound of S3 while the second line gives 
the upper bound. For y > 2, the upper and lower bounds are reversed. These two lines 
always cross with each other at a particular point, (s, S3) ~ (—4.01545, 0.062), independent 
of y. At this point, the ratio between two separations diverge (i.e., ^) and it is not 
possible to satisfy both conditions, (6.51) and (6.52), simultaneously for —4.01545 < s < 0. 
Furthermore, it can be easily checked that the upper bound of S3 is always smaller than 
the maximum value, —2 + v5, which is reached at y = and s — )■ —00. Therefore, we can 
see the physical region is confined in a very narrow semi-infinite band^^ in (s, S3)-space. 



^*In fact, even if we allow si and S2 to be different, their difference is still quite limited (i.e., |si — S2I ^ 1) 
no matter what the choice of p' is. Therefore, in general, the physical region is confined in a very narrow 
semi-inifite box. 
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Also, we carefully checked that the strong existence conditions are indeed satisfied in the 
physical region by using a similar procedure to the previous subsection. 

To understand the physical interpretation for this sector, we verify the ratio, z^, in 
(6.49). Interestingly, unlike sector (I), the ratio in the physical region is actually bounded 
above for < y < 2: 

y = l(-l-4ss3 + si) < \{-l+'^). (6.55) 

In order to have a large ratio, y must be close to zero. Since y controls how large the 
negative p^ is compared to the remaining p , this implies this negative p^ must be very 
small compared to the sum of the three p^ in order to have the solution resemble the 
classical zero entropy black ring. On the other hand, for y > 2, the maximum value of y is 
^ and moreover, the sum of the three p^ becomes negative. In this case, even though the 
solutions exist, they are no longer similar to classical black rings. 

To understand the physical interpretation of the solutions when y is not so small, 
we investigate another possibility. Specifically, we find that three centers are very likely to 
form scaling solutions in this sector. To study this possibility, we investigate the conditions 
for the formation of the triangular scaling solutions: 

(Fa , Ff,) (Ffe , Fq) > , and two other cyclic ones , (6.56) 

I (Fa , Ff,) + (F;, , Fq) I > | (Fq , Fq) | , and another two cyclic ones , (6.57) 

where the first three conditions require the directions of the fluxes to form a closed loop 
and the next three conditions require them to satisfy triangle inequalities. These conditions 
can be reduced to ones that involve only s and S3. If we look at the region, s < — 1 and 
S3 > 0, the triangle scaling conditions are reduced to: 



-5 + So -4 - S3 + s^ 

4S3 < ' < 2(l + s2 



< s < on , 2N > 0<S3<1. (6.58) 

si 



Unlike the physical region of the split flow component, this triangle scaling region defined 
by the above conditions does not depend on the dipole charges. Therefore, for any positive 
y, the scaling region remains the same. On the other hand, the physical region changes 
with y and we compare the two regions for several values of y in Figure 16. 

We find that for y > 0.4, the physical region is completely submerged in the triangle 
scaling region. For a larger value of y, the scaling region dominates and this implies that 
the three centers are more likely to form a scaling solution than become a three-center 
tree. Additionally, we randomly pick several solutions in the region where the physical and 
scaling regions overlap, we find that those solutions have bifurcation points similar to the 
example in section 5.4. This means that the split flow component only appears from the 
scaling component at the bifurcation point which is usually close to the split point. If the 
solution becomes scaling and is well-behaved, it is in fact more similar to a flve dimensional 
black hole rather than a black ring. Since a triangular scaling solution is continuously 
connected to a scaling point, based on the dynamic version of the conjecture, we can check 
Q and V Zj along the single center flow of the total charge to ensure it is well-behaved. 
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Figure 16. These graphics compare two regions in (s, s^^j-space for four values of y. The light blue 
region is the physical region of the tree while the light green one is the region that satisfies triangle 
scaling conditions. 

We find that as long as the asymptotic charges are positive (i.e., s < — 1 and S3 > —1), 
every coefficient of the radial expansion of every CTC's condition is positive. Therefore, 
the solutions in the triangle scaling region are well-behaved and should be regarded as the 
microstate geometries of some five dimensional black holes. 

Prom the above analysis, we conclude that the majority of sector (II), especially for 
the solutions with larger y, should be regarded as the microstate geometries of some black 
holes while the solutions with smaller y can be regarded as the transition phase from black 
rings to black holes. 

6.3 Sector (III) 

For this sector, we cannot find an efficient way to parametrize the charge parameters. 
However, through some numerical studies, we know the allowed region in the charge space is 



also very limited. For example, if one takes equal charge case (i.e., fc^ ^ = ^^ fe 



fc?, 



nin 



.,b ~ '^a,b ~ '^a,b ~ '^a,b), 

one finds that no tree can exist. However, one can still find some examples in which the 
tree does exist. A particular simple example is: 



-.(m) 



(-1,0,0,0), ri'") = ( 1, 150, 40, 20) , rp = (l, 10, 50, 30 



(6.59) 



We find that the tree exists in this setting and the corresponding solution are well-behaved. 
The asymptotic charges for this solution are all positive but the physical interpretation is 
unclear. 
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7 Conclusions 

In this paper, we have proposed a procedure to clarify part of CTC's free sector in bubble 
geometries based on ideas similar to the split attractor flow conjecture. Using this proce- 
dure and starting from any particular set of centers including the asymptotic vector, one 
can derive every existing tree without referring to the moduli space. The structure of the 
tree is determined by the simple thin graph that we called a "skeleton" and the existence 
conditions are based purely on checking CTC's conditions at several points or lines on this 
skeleton. If the conjecture is true, every existing tree should correspond to a split flow 
component in which every solution is well-behaved. 

By studying several numerical examples, we have shown how this correspondence ex- 
plicitly works through two pictures with one involving locating the "image" of the flows in 
space and the other involving disassembling and assembling of the solutions. We also find 
that the conjecture works well for the split flow solutions and we cannot find any counter 
example. Although, this procedure seems to be quite effective even with some complicated 
trees, it is also clear that this is not the full story. The scaling components in which some or 
all of the centers form the scaling solutions cannot be captured by this procedure. Through 
the numerical studies, we show explicitly that the scaling solutions behave like some single- 
center solutions in the split fiow pictures even though they are composed of several GH 
centers. On the other hand, in some numerical examples, we also find that the ability to 
distinguish between the scaling and split flow components may depend on the asymptotic 
vector (moduli). Specifically, when the split flow and scaling components coexist, the two 
components may exist as a whole component and only split into two when the asymptotic 
vector is tuned to the bifurcation point. It is still unclear if this is a general behavior or 
just some exotic phenomena. Also, this bifurcation point may have some physical meaning 
which is unknown at this moment. 

We also associate the simplest bubbled supertubes with some three-center trees and 
see how the existence conditions of the trees constrain their charge parameter space. The 
resulting region yields the CTC's free sector of the simplest bubbled supertubes based on 
the conjecture. Inside this region, the majority of the solutions indeed look like classical 
zero-entropy black rings. However, for solutions near the boundary of the region, the flux 
that supports the radius of the ring becomes too small and makes the solutions become 
more like black holes rather than rings. In fact, near the boundary or even some region 
beyond it, the charge parameter space becomes dominated by the scaling region and the 
three centers are more likely to form a scaling solution than a split flow solution. If they 
do form the scaling solution, this system has the total GH charge equal to one and should 
be considered as microstate geometries of some black holes according to [8] . There exists 
an overlap region for which the physical interpretation is ambiguous. A simple possibility 
is to consider this region as the transition phase between black holes and rings. 

Through several numerical examples, it is clear that the scaling solution plays an 
important role in microstate geometries of large black rings/holes. Speciflcally, it seems 
that the only way to produce black-hole like attractor points in bubble geometries is through 
the scaling solutions. However, the procedure in this paper cannot capture these solutions 
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and a similar systematic method to analyzing them is needed. If the bifurcation point we 
have discovered is quite general, it may provide a hint to search the scaling solutions by 
using the split flow trees. 

Furthermore, it is also interesting to understand what exactly happens in five dimen- 
sional moduli space. In [22], it was found that the attractors of the black rings and the 
black holes are determined from the extremization of the two different central charge func- 
tions, Zra and Zg. These functions are governed by the dipole charges and asymptotic 
electric charges respectively. Furthermore, for the attractor flow of a black ring, neither 
function is monotonic. While Zf. decreases for a large radius, it is the function Z^ which 
decreases for a small radius as one approaches the ring from infinity. This strange behavior 
is actually quite natural from the perspective of the split flow picture. The attractor flows 
simply "switch" from the flow with GH charge one to the flow with GH charge zero after 
passing the split point and the crossover region is naturally near the split point. From this 
picture, one may be able to find a unique expression for the attractor point and central 
charge for a single flow and understand the corresponding split attractor flow picture in 
five dimensional moduli space. We leave this issue and the problem related with the scaling 
solutions to a future study. 
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